rownames(res98_full_tab) <- c("Intercept", "Treatment", "CA PCA score", "Treatment X CA PCA score",
"HAND_kapisa", "HAND_parwan", "cv_ts")
tab_a12_ts_ca_pca_full<-multimerge( list (one=res95_full_tab, two=res96_full_tab, three=res97_full_tab, four=res98_full_tab))
tab_a12_ts_ca_pca_full_format_print<-xtable(tab_a12_ts_ca_pca_full,
caption="\\label{tab:tab_a12_ts_ca_pca_full}Interaction Effects Regressions with Village-Level PCA Scores of Current CDC Functionality with Agency Outcomes for Teachers")
columns<-c("TS Fund Application Index", "TS Classroom Supplies Score", "TS Motivation Index", "TS Paid Regularly")
N<-c(length(lm95$residuals), length(lm96$residuals), length(lm97$residuals), length(lm98$residuals))
addtorow<- list()
addtorow$pos <-list(0)
addtorow$command <-paste0(paste0('& \\multicolumn{2}{c}{', unique(columns), '}', collapse=''), '\\\\',
paste0('&\\multicolumn{2}{c}{', N, '}', collapse=''), '\\\\')
sink(file="tabs-figs/tab_a12_ts_ca_pca_full.tex")
print(tab_a12_ts_ca_pca_full_format_print,
caption.placement = "top",
add.to.row=addtorow,
hline.after = -1)
sink()
##########################################################
#### FULL COVARIATE OUTPUT FOR TABLE A13 IN MAIN APPENDIX
##########################################################
# one table for the community leader outcomes
#combined sustainability index
res30_full_tab<- data.frame(Coef. = round(res30_full$beta, 2),
SE = sapply(round(res30_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res30_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#community provision for future classes
res31_full_tab<- data.frame(Treat.Effect  = round(res31_full$beta, 2),
SE = sapply(round(res31_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res31_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#confidence in community institutions
res32_full_tab<- data.frame(Treat.Effect  = round(res32_full$beta, 2),
SE = sapply(round(res32_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res32_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#confidence that school will be sustained
res33_full_tab<- data.frame(Treat.Effect  = round(res33_full$beta, 2),
SE = sapply(round(res33_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res33_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#MOE provisions for funding teachers
res34_full_tab<- data.frame(Treat.Effect  = round(res34_full$beta, 2),
SE = sapply(round(res34_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res34_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Shura provisions for future classes
res35_full_tab<- data.frame(Treat.Effect  = round(res35_full$beta, 2),
SE = sapply(round(res35_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res35_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
tab_a13_cl_full<-multimerge( list (one=res30_full_tab, two=res31_full_tab, three=res32_full_tab,
four=res33_full_tab, five=res34_full_tab, six=res35_full_tab))
tab_a13_cl_full_format_print<-xtable(tab_a13_cl_full,
caption="\\label{tab:tab_a13_cl_full}Effects of Sustainability Model on Sustainability Outcomes for Community Leaders")
columns<-c("Combined Sustainability Index", "Community Provision for Future Classes", "Confidence in Community Instns.",
"Confidence that School Will be Sustained", "MOE Provisions for Funding Teachers", "Shura Provisions for Future Classes")
N<-c(length(lm30$residuals), length(lm31$residuals), length(lm32$residuals),
length(lm33$residuals), length(lm34$residuals), length(lm35$residuals))
addtorow<- list()
addtorow$pos <-list(0)
addtorow$command <-paste0(paste0('& \\multicolumn{2}{c}{', unique(columns), '}', collapse=''), '\\\\',
paste0('&\\multicolumn{2}{c}{', N, '}', collapse=''), '\\\\')
sink(file="tabs-figs/tab_a13_cl_full.tex")
print(tab_a13_cl_full_format_print,
caption.placement = "top",
add.to.row=addtorow,
hline.after = -1)
sink()
#includes a little manual formatting
# one table for the teacher outcomes
#Combined sustainability index index
res36_full_tab<- data.frame(Coef. = round(res36_full$beta, 2),
SE = sapply(round(res36_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res36_full_tab) <- c("Intercept", "Treatment", "ts_cv")
#confidence in local community institutions
res37_full_tab<- data.frame(Coef.  = round(res37_full$beta, 2),
SE = sapply(round(res37_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res37_full_tab) <- c("Intercept", "Treatment", "ts_cv")
#Confidence that school will be sustained
res38_full_tab<- data.frame(Coef. = round(res38_full$beta, 2),
SE = sapply(round(res38_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res38_full_tab) <- c("Intercept", "Treatment", "ts_cv")
#proactive management from MOE
res39_full_tab<- data.frame(Coef. = round(res39_full$beta, 2),
SE = sapply(round(res39_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res39_full_tab) <- c("Intercept", "Treatment", "ts_cv")
#Proactive management from shuras
res40_full_tab<- data.frame(Coef. = round(res40_full$beta, 2),
SE = sapply(round(res40_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res40_full_tab) <- c("Intercept", "Treatment", "ts_cv")
#Teacher plans to remain CBE teacher
res41_full_tab<- data.frame(Coef.  = round(res41_full$beta, 2),
SE = sapply(round(res41_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res41_full_tab) <- c("Intercept", "Treatment", "ts_cv")
tab_a13_ts_full<-multimerge( list (one=res36_full_tab, two=res37_full_tab, three=res38_full_tab,
four=res39_full_tab, five=res40_full_tab, six=res41_full_tab))
tab_a13_ts_full_format_print<-xtable(tab_a13_ts_full,
caption="\\label{tab:tab_a13_ts_full}Effects of Sustainability Model on Sustainability Outcomes for Teachers")
columns<-c("Combined Sustainability Index", "Confidence in Local Community Instns.", "Confidence that School Will be Sustained",
"Proactive Management from MOE", "Proactive Management from Shuras", "Teacher plans to remain CBE Teacher")
N<-c(length(lm36$residuals), length(lm37$residuals), length(lm38$residuals),
length(lm39$residuals), length(lm40$residuals), length(lm41$residuals))
addtorow<- list()
addtorow$pos <-list(0)
addtorow$command <-paste0(paste0('& \\multicolumn{2}{c}{', unique(columns), '}', collapse=''), '\\\\',
paste0('&\\multicolumn{2}{c}{', N, '}', collapse=''), '\\\\')
sink(file="tabs-figs/tab_a13_ts_full.tex")
print(tab_a13_ts_full_format_print,
caption.placement = "top",
add.to.row=addtorow,
hline.after = -1)
sink()
#includes a little manual formatting
##########################################################
#### FULL COVARIATE OUTPUT FOR TABLE A14 IN MAIN APPENDIX
##########################################################
#  child outcomes and covariates
# attendance
res1.crs_full_tab<- data.frame(Coef. = round(res1.crs_full$beta, 2),
SE = sapply(round(res1.crs_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res1.crs_full_tab) <- c("Intercept", "Treatment", "CRS", "Treatment X CRS",
"ros_hhheadchild", "ros_girl", "ros_age", "int_lang_pashto", "hh_eth_hazara", "hh_eth_pashtun", "hh_eth_tajik",
"hhhead_job_farmer", "hhhead_age", "hhhead_schoolyears", "hh_num", "hh_jeribs", "hh_animal_sheep", "school_km",
"hhhead_job_laborer", "hhhead_job_military", "hhhead_write", "hhhead_read", "hhhead_edu_mosque", "hhhead_edu_madrassa",
"hhhead_edu_community", "hhhead_edu_government", "hhhead_edu_university", "hhhead_female", "hh_totalincome_2001less",
"hh_totalincome_2001to5000", "hhhead_totalincome_5001to10000", "hhhead_totalincome_10001to15000", "hh_totalincome_15001plus",
"hh_own_tvs", "hh_own_mobiles", "hh_own_cars", "hh_own_radios", "factorwealth", "hh_landown", "hh_children", "PROVINCE1",
"PROVINCE3", "PROVINCE4", "PROVINCE5")
# test scores
res2.crs_full_tab<- data.frame(Coef. = round(res2.crs_full$beta, 2),
SE = sapply(round(res2.crs_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res2.crs_full_tab) <- c("Intercept", "Treatment", "CRS", "Treatment X CRS",
"ros_hhheadchild", "ros_girl", "ros_age", "int_lang_pashto", "hh_eth_hazara", "hh_eth_pashtun", "hh_eth_tajik",
"hhhead_job_farmer", "hhhead_age", "hhhead_schoolyears", "hh_num", "hh_jeribs", "hh_animal_sheep", "school_km",
"hhhead_job_laborer", "hhhead_job_military", "hhhead_write", "hhhead_read", "hhhead_edu_mosque", "hhhead_edu_madrassa",
"hhhead_edu_community", "hhhead_edu_government", "hhhead_edu_university", "hhhead_female", "hh_totalincome_2001less",
"hh_totalincome_2001to5000", "hhhead_totalincome_5001to10000", "hhhead_totalincome_10001to15000", "hh_totalincome_15001plus",
"hh_own_tvs", "hh_own_mobiles", "hh_own_cars", "hh_own_radios", "factorwealth", "hh_landown", "hh_children", "PROVINCE1",
"PROVINCE3", "PROVINCE4", "PROVINCE5")
#  household satisfaction outcome and covariates
# test scores
res3.crs_full_tab<- data.frame(Coef. = round(res3.crs_full$beta, 2),
SE = sapply(round(res3.crs_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res3.crs_full_tab) <- c("Intercept", "Treatment", "CRS", "Treatment X CRS",
"int_lang_pashto", "hh_eth_hazara", "hh_eth_pashtun", "hh_eth_tajik",
"hhhead_job_farmer", "hhhead_age", "hhhead_schoolyears", "hh_num", "hh_jeribs", "hh_animal_sheep", "school_km",
"hhhead_job_laborer", "hhhead_job_military", "hhhead_write", "hhhead_read", "hhhead_edu_mosque", "hhhead_edu_madrassa",
"hhhead_edu_community", "hhhead_edu_government", "hhhead_edu_university", "hhhead_female", "hh_totalincome_2001less",
"hh_totalincome_2001to5000", "hhhead_totalincome_5001to10000", "hhhead_totalincome_10001to15000", "hh_totalincome_15001plus",
"hh_own_tvs", "hh_own_mobiles", "hh_own_cars", "hh_own_radios", "factorwealth", "hh_landown", "hh_children", "PROVINCE1",
"PROVINCE3", "PROVINCE4", "PROVINCE5")
tab_a14_full<-multimerge( list (one=res1.crs_full_tab, two=res2.crs_full_tab, three=res3.crs_full_tab))
tab_a14_full_format_print<-xtable(tab_a14_full,
caption="\\label{tab:tab_a14_full}Heterogeneous Treatment Effects by Implementing NGO")
columns<-c("Child Attendance", "Child Test Score", "Household Satisfaction")
N<-c(length(lm1.crs$residuals), length(lm2.crs$residuals), length(lm3.crs$residuals))
addtorow<- list()
addtorow$pos <-list(0)
addtorow$command <-paste0(paste0('& \\multicolumn{2}{c}{', unique(columns), '}', collapse=''), '\\\\',
paste0('&\\multicolumn{2}{c}{', N, '}', collapse=''), '\\\\')
sink(file="tabs-figs/tab_a14_full.tex")
print(tab_a14_full_format_print,
caption.placement = "top",
add.to.row=addtorow,
hline.after = -1)
sink()
#requires a little bit of manual formatting, added in het.p-value from main table
##########################################################
#### FULL COVARIATE OUTPUT FOR TABLE A15 IN MAIN APPENDIX
##########################################################
#one table for the equivalency outcome
res42_full_tab<- data.frame(Coef.= round(res42_full$beta, 2),
SE = sapply(round(res42_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res42_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
res42_full_tab_format_print<-xtable(res42_full_tab,
caption="\\label{tab:tab_a15_full_equiv}Community Leaders' Longer-Term Perceptions of CBE Sustainability (CARE and CRS Samples) - Equivalency Outcome")
columns<-c("Community Desire to Continue CBE")
N<-c(length(lm42$residuals))
addtorow<- list()
addtorow$pos <-list(0)
addtorow$command <-paste0(paste0('& \\multicolumn{2}{c}{', unique(columns), '}', collapse=''), '\\\\',
paste0('&\\multicolumn{2}{c}{', N, '}', collapse=''), '\\\\')
sink(file="tabs-figs/tab_a15_full_equiv.tex")
print(res42_full_tab_format_print,
caption.placement = "top",
add.to.row=addtorow,
hline.after = -1)
sink()
#one table for the non-equivalency outcomes
# combined sustainability outcome
res43_full_tab<- data.frame(Coef. = round(res43_full$beta, 2),
SE = sapply(round(res43_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res43_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Initiative to have CBE sustained
res44_full_tab<- data.frame(Coef. = round(res44_full$beta, 2),
SE = sapply(round(res44_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res44_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Provisions for sustainability without NGO
res45_full_tab<- data.frame(Coef. = round(res45_full$beta, 2),
SE = sapply(round(res45_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res45_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Confidence that school will sustain
res46_full_tab<- data.frame(Coef. = round(res46_full$beta, 2),
SE = sapply(round(res46_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res46_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Continued CBE operations in spring 2018
res47_full_tab<- data.frame(Coef. = round(res47_full$beta, 2),
SE = sapply(round(res47_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res47_full_tab) <- c("Intercept", "Treatment", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
tab_a15_full_nonequiv<-multimerge( list (one=res43_full_tab, two=res44_full_tab, three=res45_full_tab, four=res46_full_tab, five=res47_full_tab))
tab_a15_full_nonequiv_format_print<-xtable(tab_a15_full_nonequiv,
caption="\\label{tab:tab_a15_full_nonequiv}Community Leaders' Longer-Term Perceptions of CBE Sustainability (CARE and CRS Samples) - Non-Equivalency Outcomes")
columns<-c("Combined Sustainability Index", "Initiative to have CBE Sustained", "Provisions for Sustainability Without NGO",
"Confidence that School Will Sustain", "Confinued CBE Operations in Spring 2018")
N<-c(length(lm44$residuals), length(lm44$residuals), length(lm45$residuals), length(lm46$residuals), length(lm47$residuals))
addtorow<- list()
addtorow$pos <-list(0)
addtorow$command <-paste0(paste0('& \\multicolumn{2}{c}{', unique(columns), '}', collapse=''), '\\\\',
paste0('&\\multicolumn{2}{c}{', N, '}', collapse=''), '\\\\')
sink(file="tabs-figs/tab_a15_full_nonequiv.tex")
print(tab_a15_full_nonequiv_format_print,
caption.placement = "top",
add.to.row=addtorow,
hline.after = -1)
sink()
#need to do a little manul formatting
##########################################################
#### FULL COVARIATE OUTPUT FOR TABLE A16 IN MAIN APPENDIX
##########################################################
#one table with the regressions for community leader's age
#combined sustainability index
res48_full_tab<- data.frame(Coef.= round(res48_full$beta, 2),
SE = sapply(round(res48_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res48_full_tab) <- c("Intercept", "Treatment", "Older leader",
"Treatment x Older Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Community provision for future classes
res48b_full_tab<- data.frame(Coef.= round(res48b_full$beta, 2),
SE = sapply(round(res48b_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res48b_full_tab) <- c("Intercept", "Treatment", "Older leader",
"Treatment x Older Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Confidence in community institutions
res48c_full_tab<- data.frame(Coef.= round(res48c_full$beta, 2),
SE = sapply(round(res48c_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res48c_full_tab) <- c("Intercept", "Treatment", "Older leader",
"Treatment x Older Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Confidence that school will be sustained
res48d_full_tab<- data.frame(Coef.= round(res48d_full$beta, 2),
SE = sapply(round(res48d_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res48d_full_tab) <- c("Intercept", "Treatment", "Older leader",
"Treatment x Older Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#MOE provisions for funding teachers
res48e_full_tab<- data.frame(Coef.= round(res48e_full$beta, 2),
SE = sapply(round(res48e_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res48e_full_tab) <- c("Intercept", "Treatment", "Older leader",
"Treatment x Older Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#MOE provisions for funding teachers
res48f_full_tab<- data.frame(Coef.= round(res48f_full$beta, 2),
SE = sapply(round(res48f_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res48f_full_tab) <- c("Intercept", "Treatment", "Older leader",
"Treatment x Older Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
tab_a16_full_age<-multimerge( list (one=res48_full_tab, two=res48b_full_tab, three=res48c_full_tab,
four=res48d_full_tab, five=res48e_full_tab, six=res48f_full_tab))
tab_a16_full_age_format_print<-xtable(tab_a16_full_age,
caption="\\label{tab:tab_a16_full_age}Interaction Effect Regressions With Community Leaders' Age, Endline Community Leader Sustainability Perceptions")
columns<-c("Combined Sustainability Index", "Community Provisions for Future Classes", "Confidence in Community Instns.",
"Confidence that School Will Be Sustained", "MOE provisions for Funding Teachers", "Shura Provisions for Future Classes")
N<-c(length(lm48$residuals), length(lm48b$residuals), length(lm48c$residuals), length(lm48d$residuals), length(lm48e$residuals), length(lm48f$residuals))
addtorow<- list()
addtorow$pos <-list(0)
addtorow$command <-paste0(paste0('& \\multicolumn{2}{c}{', unique(columns), '}', collapse=''), '\\\\',
paste0('&\\multicolumn{2}{c}{', N, '}', collapse=''), '\\\\')
sink(file="tabs-figs/tab_a16_full_age.tex")
print(tab_a16_full_age_format_print,
caption.placement = "top",
add.to.row=addtorow,
hline.after = -1)
sink()
#one table with the regressions for community leader's education
#combined sustainability index
res49_full_tab<- data.frame(Coef.= round(res49_full$beta, 2),
SE = sapply(round(res49_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res49_full_tab) <- c("Intercept", "Treatment", "Educated leader",
"Treatment x Educated Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Community provision for future classes
res49b_full_tab<- data.frame(Coef.= round(res49b_full$beta, 2),
SE = sapply(round(res49b_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res49b_full_tab) <- c("Intercept", "Treatment", "Educated leader",
"Treatment x Educated Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Confidence in community institutions
res49c_full_tab<- data.frame(Coef.= round(res49c_full$beta, 2),
SE = sapply(round(res49c_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res49c_full_tab) <- c("Intercept", "Treatment", "Educated leader",
"Treatment x Educated Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Confidence that school will be sustained
res49d_full_tab<- data.frame(Coef.= round(res49d_full$beta, 2),
SE = sapply(round(res49d_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res49d_full_tab) <- c("Intercept", "Treatment", "Educated leader",
"Treatment x Educated Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#MOE provisions for funding teachers
res49e_full_tab<- data.frame(Coef.= round(res48e_full$beta, 2),
SE = sapply(round(res48e_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res49e_full_tab) <- c("Intercept", "Treatment", "Educated leader",
"Treatment x Educated Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#MOE provisions for funding teachers
res49f_full_tab<- data.frame(Coef.= round(res49f_full$beta, 2),
SE = sapply(round(res49f_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res49f_full_tab) <- c("Intercept", "Treatment", "Educated leader",
"Treatment x Educated Leader", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
tab_a16_full_educ<-multimerge( list (one=res49_full_tab, two=res49b_full_tab, three=res49c_full_tab,
four=res49d_full_tab, five=res49e_full_tab, six=res49f_full_tab))
tab_a16_full_educ_format_print<-xtable(tab_a16_full_educ,
caption="\\label{tab:tab_a16_full_educ}Interaction Effect Regressions With Community Leaders' Educational Status, Endline Community Leader Sustainability Perceptions")
columns<-c("Combined Sustainability Index", "Community Provisions for Future Classes", "Confidence in Community Instns.",
"Confidence that School Will Be Sustained", "MOE provisions for Funding Teachers", "Shura Provisions for Future Classes")
N<-c(length(lm49$residuals), length(lm49b$residuals), length(lm49c$residuals), length(lm49d$residuals), length(lm49e$residuals), length(lm49f$residuals))
addtorow<- list()
addtorow$pos <-list(0)
addtorow$command <-paste0(paste0('& \\multicolumn{2}{c}{', unique(columns), '}', collapse=''), '\\\\',
paste0('&\\multicolumn{2}{c}{', N, '}', collapse=''), '\\\\')
sink(file="tabs-figs/tab_a16_full_educ.tex")
print(tab_a16_full_educ_format_print,
caption.placement = "top",
add.to.row=addtorow,
hline.after = -1)
sink()
##########################################################
#### FULL COVARIATE OUTPUT FOR TABLE A17 IN MAIN APPENDIX
##########################################################
#combined sustainability index
res106_full_tab<- data.frame(Coef.= round(res106_full$beta, 2),
SE = sapply(round(res106_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res106_full_tab) <- c("Intercept", "Treatment", "Mean Test Score",
"Treatment x Mean Test Score", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Community provision for future classes
res101_full_tab<- data.frame(Coef.= round(res101_full$beta, 2),
SE = sapply(round(res101_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res101_full_tab) <- c("Intercept", "Treatment", "Mean Test Score",
"Treatment x Mean Test Score", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Confidence in community institutions
res102_full_tab<- data.frame(Coef.= round(res102_full$beta, 2),
SE = sapply(round(res102_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res102_full_tab) <- c("Intercept", "Treatment", "Mean Test Score",
"Treatment x Mean Test Score", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Confidence that school will be sustained
res103_full_tab<- data.frame(Coef.= round(res103_full$beta, 2),
SE = sapply(round(res103_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res103_full_tab) <- c("Intercept", "Treatment", "Mean Test Score",
"Treatment x Mean Test Score", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#MOE provisions for funding teachers
res104_full_tab<- data.frame(Coef.= round(res104_full$beta, 2),
SE = sapply(round(res104_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res104_full_tab) <- c("Intercept", "Treatment", "Mean Test Score",
"Treatment x Mean Test Score", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
#Shura provisions for future classes
res105_full_tab<- data.frame(Coef.= round(res105_full$beta, 2),
SE = sapply(round(res105_full$SE, 2),
function(x) {paste0("(",x,")")}))
rownames(res105_full_tab) <- c("Intercept", "Treatment", "Mean Test Score",
"Treatment x Mean Test Score", "int_lang_pashto", "hhhead_job_farmer",
"hhhead_age", "hhhead_schoolyears", "school_km")
tab_a17_full<-multimerge( list (one=res106_full_tab, two=res101_full_tab, three=res102_full_tab,
four=res103_full_tab, five=res104_full_tab, six=res105_full_tab))
tab_a17_full_format_print<-xtable(tab_a17_full,
caption="\\label{tab:tab_a17_full}Interaction Effect Regressions With Community Mean Children's Test Score, Endline Community Leader Sustainability Perceptions")
columns<-c("Combined Sustainability Index", "Community Provisions for Future Classes", "Confidence in Community Instns.",
"Confidence that School Will Be Sustained", "MOE provisions for Funding Teachers", "Shura Provisions for Future Classes")
N<-c(length(lm106$residuals), length(lm101$residuals), length(lm102$residuals), length(lm103$residuals), length(lm104$residuals), length(lm105$residuals))
addtorow<- list()
addtorow$pos <-list(0)
addtorow$command <-paste0(paste0('& \\multicolumn{2}{c}{', unique(columns), '}', collapse=''), '\\\\',
paste0('&\\multicolumn{2}{c}{', N, '}', collapse=''), '\\\\')
sink(file="tabs-figs/tab_a17_full.tex")
print(tab_a17_full_format_print,
caption.placement = "top",
add.to.row=addtorow,
hline.after = -1)
sink()
#############################################################################################
### TABLE A22: EFFECT OF TREATMENT ON PARENTAL INVOLVEMENT
############################################################################################
## Household outcome 1: parental involvement in student's homework
y.p1 <- within_hh_dat$q37_recode
lmp1<-lm( formula(paste("y.p1", "XX_hh", sep="~")),  weights=within_hh_dat$wt)
resp1<-coef_test(lmp1, vcov=vcovCR(lmp1, cluster=within_hh_dat$CLUSTER, type="CR2"))
print(resp1)
#full results with all covariates [for appendix]
resp1full<-as.data.frame(resp1)
#main results without all the covariates [for main body of paper]
resp1<-as.data.frame(resp1[2,]) # keep second row of output
# add in columns on the sesoi for teacher_attend_score
sesoip1<- sesoi_gen
zp1<-(resp1[,"beta"]-sesoip1)/resp1[,"SE"]
p_sesoip1=pt(-zp1, df = resp1[,"df_Satt"])
# add in t_Satt_citical
tp1 <-sapply(resp1$df, function(d) qt(1-(1-0.9)/2, df=d)) # ci= 0.9
# final set of results for parental involvement in homework
resp1<-cbind(resp1, sesoi=sesoip1, p_sesoi=p_sesoip1, t_Satt_crit=tp1, ci_l=resp1$beta-tp1*resp1[,"SE"],
ci_u=resp1[,"beta"]+tp1*resp1[,"SE"],
coverage=0.9, N=length(lmp1$residuals))
rownames(resp1)<-paste(rownames(resp1), 'parent_involv_hw')
resp1
# Household outcome 2: parental contributions
y.p2<-within_hh_dat$q40
lmp2<-lm( formula(paste("y.p2", "XX_hh", sep="~")),  weights=within_hh_dat$wt)
resp2<-coef_test(lmp2, vcov=vcovCR(lmp2, cluster=within_hh_dat$CLUSTER, type="CR2"))
print(resp2)
#full results with all covariates [for appendix]
resp2full<-as.data.frame(resp2)
#main results without all the covariates [for main body of paper]
resp2<-as.data.frame(resp2[2,]) # keep second row of output
# add in columns on the sesoi for parent_contribution
sesoip2<- sesoi_gen
zp2<-(resp2[,"beta"]-sesoip2)/resp2[,"SE"]
p_sesoip2=pt(-zp2, df = resp2[,"df_Satt"])
# add in t_Satt_citical
tp2 <-sapply(resp2$df, function(d) qt(1-(1-0.9)/2, df=d)) # ci= 0.9
# final set of results for parent_contribution
resp2<-cbind(resp2, sesoi=sesoip2, p_sesoi=p_sesoip2, t_Satt_crit=tp2, ci_l=resp2$beta-tp2*resp2[,"SE"],
ci_u=resp2[,"beta"]+tp2*resp2[,"SE"],
coverage=0.9, N=length(lmp2$residuals))
rownames(resp2)<-paste(rownames(resp2), 'parental_contributions')
resp2
# combining results in Table together
statp1<-jmsummary(within_hh_dat$q37_recode)
rownames(statp1) <- "Parental involvement in homework"
statp2<-jmsummary(within_hh_dat$q40)
rownames(statp2) <- "Parental contributions"
TableX_Output<-rbind(resp1, resp2)
TableX_sum <- rbind(statp1, statp2)
TableX_format <- data.frame(Treatment.Effect = round(TableX_Output$beta, 2),
Std.Err. = sapply(round(TableX_Output$SE, 2),
function(x) {paste0("(",x,")")}),
N = TableX_Output$N,
Control.Mean = round(TableX_sum[,"mean"],2),
Control.Std.Dev. = round(TableX_sum[,"sd"],2))
rownames(TableX_format) <- rownames(TableX_sum)
TableX_print<-xtable(TableX_format,
caption="\\label{tab:parents_involv}Parent's Involvement and Contributions")
addtorow <- list()
addtorow$pos <- list(0,1,1)
addtorow$command <- c("{\  &  &  &  & &  & &\\\\\n",
"\\\\\n",
"{\  &  &  &  & &  & & \\\\\n")
sink(file="tabs-figs/hh_parent_involvment.tex")
print(TableX_print,
add.to.row = addtorow,
caption.placement = "top",
hline.after = -1)
sink()
table(within_hh_dat$q40)
